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Abstract. We consider the hierarchy of higher-order Riccati equations and establish their 
connection with the Gambier equation. Moreover we investigate the relation of equations 
of the Gambier family to other nonlinear differential systems. In particular we explore 
their connection to the generalized Ermakov-Pinney and Milne-Pinncy equations. In addi- 
tion we investigate the consequence of introducing Okamoto's folding transformation which 
maps the reduced Gambier equation to a Lienard type equation. Finally the conjugate 
Hamiltonian aspects of certain equations belonging to this family and their connection with 
superintegrability are explored. 
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1 Introduction 



The Gambier equation plays, for linearizable systems, the same role as the Painleve VI. Just 
as one can derive, starting from the latter, all the remaining Painleve equations through appro- 
priate limiting procedures, in the same way the Gambier equation is a "master" system for the 
linearizable equations of the Painleve-Gambier classification \17\ [HI [13] . The linearizability of 
the Gambier equation is based on the fact that it consists in two Riccati equations in cascade 
(where by "cascade" we mean that the solution of the first Riccati equation is used in order to 
compute the coefficients of the second one). Thus the Gambier equation, written as a system, 
has the form: 



-y +by + c, 
ax^ + nyx -\- a. 



(1.1a) 
(1.1b) 



Here a, h and c are functions of the independent variable, o" is a constant which can be scaled 
to 1 unless it happens to be and n is an integer. Eliminating y between (jl.lap and (jl.lbp one 



*This paper is a contribution to the Proceedings of the Conference "Integrable Systems and Geomet- 
ry" (August 12-17, 2010, Pondicherry University, Puducherry, India). The full collection is available at 
ihttp://www.emis.de/journals/SIGMA/ISG2010.html| 
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can write the Gambier equation [TOl [TT] as a second-order equation: 

n — 1 iP' n + 2 n — 2 x a? n 

X = h a XX + hx (T X 

n X n n X n 

+ (a- ah)x^ + (cn-'^^\x -ba - —. (1.2) 
V n J nx 



Equation (|1.2p . although the most general one, is not the only one in the Gambier list. We have 
in fact 

X = —3xx — + qi^x + x^) , (G5) 

.2, (G13) 

,x3 + ^x^ (G14) 











X 




— + 


q- 


- q 


x 


X 






i) 


x? 


n 


X 



(n + 2)2 n + 2 

where we have preserved the numbering of Gambier [8]. To this Ust of Unearizable systems one 
must in principle adjoin the equation 

X = —2xx + qx + qx, (G6) 

which is just the derivative of the Riccati equation. 

Starting from the Gambier equation one can obtain almost all the Unearizable equations of 
the Painleve-Gambier list by the appropriate limits. We start by taking n = 1 and a = 
in (|1.2p . By the appropriate independent variable transformation and gauge we can put c = 
and a = — 1 leading to ()G5p . Next we start from (|1.2p and take the limit n — )■ oo. We find 

X 



■ 2 

X . , . X 



+ axx + bx cr + (a — ab)x'^ + dx — ba (1-3) 



(where we have introduced d = lim cn). Equation (jl.3p is not in canonical form. We can 

ri —>oo 

reduce it to the standard expression by taking 6 = and introducing a gauge x ^ px such that 
d = pI p — I p^ ■ We obtain thus precisely ()G13P where q = 1/ p and r = ap. Finally we take 
£7 = 0. and perform an independent variable transformation [11] which allows to put 6 = 0. It 
turns out that in this case an additional gauge freedom does exist allowing us to take c = 0. We 
obtain thus the equation 

n— liP'n + 2a?r,2 

X = h a XX X + ax , (1-4) 

n X n n 

which is precisely (IG14P with q = a(l + 2/n). 

The derivative of the Riccati equation, ()G6p . constitutes an exception in the sense that it 
cannot be obtained from the Gambier equation. Since we mentioned the derivative of the Riccati 
it is interesting to point out another "Riccati derivative" which is none other but (jl.Sp or, in 
canonical form, (IGlSp . Indeed at the level of the coupled Riccati for the limit n — t- oo to have 
a meaning y must go to zero as well. This means that the quadratic term in (jl.lap disappears. 
Moreover we saw that the canonical form corresponds to 6 = and since c is free, we can 
introduce a new function through d = nc. Finally, we divide (jl.lbp by x and take the derivative. 
This leads to the appearance of a term ny, which from (11. lap is equal to d. Thus, equation (11. 4p 
is a derivative of a Riccati after we have divided by the dependent variable. 

At this point one remark is in order concerning the integrability of (jl.ip . The Gambier 
equation was first derived by Gambier who complemented the study of Painleve of second order 
differential equations the solutions of which have, what is today called, the Painleve property. 
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In order to study the latter in the case of the coupled Riccati system we must start from the 
observation that since p.la|) is a Riccati the dominant behaviour of y is y ~ — zq). An 
expansion of y around this singularity can be easily established: it involves the functions b, c and 
their derivatives. Next we turn to (jl.lbp which is again a Riccati and thus its singularities are 
simple poles. However there exist singularities which are due to the singular behaviour of the 
coefficients, namely y. The locations of the singularities of the coefficients of some equation are 
considered as 'fixed', i.e. not dependent on the initial data. However, in the present case, these 
singularities are movable because the zq in y is arbitrary and thus should be studied. Because of 
the pole in y, x has a singular expansion. The expansion has a resonance at position n — 1 and 
this introduces a compatibility condition. Thus some constraint on the functions a, 6, c must 
exist for the solutions to have the Painleve property. However the latter is not necessary when 
it comes to integrating a linearisable system. The algebraic integrability which stems from the 
Painleve property may confer additional "nice" properties to the solutions of the system but, if 
what we are interested in is just to be able to integrate the system, linearizability suffices. In 
the case of the Gambier equation this means that if one contents oneself with linearizability no 
constraints on the free functions do exist and as a matter of fact n does not have to be integer. 

1.1 Motivation and result 

Curiously the Gambier equation has not been the object of many studies. In particular its 
relation to other integrable differential systems has not been investigated. This is something we 
intend to remedy in the present paper. We have several results centered around the Gambier 
equation. These are listed below. 

• We explore a connection between the generalized Ermakov-Pinney [5l [18] and the Milne- 
Pinney [Ml [18] equations with the reduced Gambier type equation. 

• Using Okamoto's folding transformation we map the reduced Gambier equation to an 
equation of the Lienard type. We also derive the conjugate Hamiltonian equations for the 
reduced Gambier equation. 

• In addition the conditions under which the reduced Gambier equation becomes a superin- 
tegrable system is investigated. 

2 Preliminaries, the generalised Riccati equation 

The properties of higher-order Riccati equations were investigated in [1] where the authors 
studied in detail the Riccati differential operator, 



and the corresponding Riccati differential equations of order m belonging to the hierarchy 



Here h is assumed to be a constant independent of t. Therefore, as a first generalization one may 
assume to be a time-dependent function and define a corresponding differential operator, Tr, 




X = 0, m = 0,1,2, ... . 




by 
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n = 


1 


Trx = 


{- 

\dt 


n = 


2 


T|x = 


i- 


n = 


3 




i- 

\dt 



given in [2]. In a manner similar to (|2.ip one may define, therefore, a nonautonomous hierarchy 
of equations by 

T^x = 0, m = 0,1,2,.... 

The first few equations of this hierarchy are given by: 

+ h{t)x^ X = X + h{t)x'^, 

+ h{t)x^ X = x + 3h{t)xx + h^{t)x^ + h(t)x'^, 

+ h{t)x^ x=x+ Ah{t)xx + Qh^{t)x^x + 3/i(t)x^ + h^{t)x'^ 

+ hh{t)xx + 3h{t)h{t)x^ + h{t)x'^, 

and analogous expressions for higher values of n. Obviously the equations become increasingly 
bulky when n > 3. Of particular interest is the second member of this family of equations, 
namely 

X + 3h{t)xx + h{t)x^ + h^{t)x^ = 0. (2.2) 

Equation (j2.2p is called the time-dependent second-order Riccati equation and is known to admit 
a Lagrangian description [2] . Recently a time-dependent first integral for this equation has been 
deduced and this has provided the impetus to investigate more deeply the nature of the Riccati 
operator Tr. 

2.1 Linearization of generalised second-order Riccati equations 
and the Ermakov— Pinney system 

It is well known that the conventional projective transformation 

X = Ij (2,3) 

with h = constant can be applied to linearize the time-independent second-order Riccati equation 
X + 3hxx + h^x^ = 0. 

When h = h{t) the same projective transformation linearizes the time-dependent second-order 
Riccati equation (12. 2p to the following, namely. 



2h ; f^fhV h 



where / denotes differentiation with respect to the argument. 

Notice that (|2.2p corresponds to the choice n = l, a = b = and a = —h hi (jl.2p . If 
6 7^ then (jl.2p reduces (after suitable redefinition) of the time dependent coefficients, namely 
b —7- —b, c —7- — c and a = —h to 

X + 3hxx + bx + cx+ [h + bh)x^ + h^x^ = 0. 

This equation has been considered by Sugai [20]. Under the projective transformation defined 
by (j2.3p its corresponding linear version is seen to be 
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2.1.1 Hill's equation and the Ermakov Pinney equation 

We wish to recap [T^] here the connection between Hill's equation and the Ermakov-Pinney 
equation, which will be generalized in the next section. It is well known that if -01 and 02 are 
the solutions of the Hill's equation 

then the product ^jJ^^pj = 1,2) satisfies the third-order linear equation 

7 + 2uf + Auf = 0, (2.4) 

and traces out a three dimensional space of solutions. Since the solution of equation (|2.4|) is 
spanned by Span(0^, 0|, 0i'(/'2)) its general solution is given by an arbitrary linear combination 
of the basis vectors, viz 

iP = AiPl + 250iV'2 + Cil^l 

One can show that this is in turn connected with 



where a is a constant. 

Proposition 1. If ipi and 11)2 satisfy Hill's equation then 



satisfies the Ermakov equation 

^ + n(t)V' = -^, cj = AC-B'^. 

Proof. Since 

= AV'? + 260^1 02 + C0i, 
differentiation yields 

+ n(t)0 = ^3 [{AiPl + 250102 + C0i) (^0? + 260102 + C0i) 

- (A0101 + 5(0102 + 0102) + C0202)^ 

Upon simplification this gives 

+ u(t)0= ^^^3^ W\ 

The proof follows by setting a = AC — and defining the Wronskian W = (0i02 — 0i02) to 
be unity. ■ 

It is possible to further extend the above result. 
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2.2 Generalized Ermakov— Pinney system and the Gambier equation 

In this section we extend the previous results. Confining ourselves to a special choice of Hill's 
equation, viz u{x) = 0^, a constant, i.e., the harmonic oscillator equation we have the following 
proposition. 

Proposition 2. Let tpi and ip2 satisfy the equation of a linear harmonic oscillator then 
satisfies 

q + n^q = -^, (2.5) 

q^ 

where K is some constant and A, B, C , D satisfy following conditions 

AC - = BD - C"^ = -A and AD - BC = 0. 

Proof. The proof proceeds in a similar manner as in Proposition [TJ Since q"^ = Atp'f + 3BipfiJj2 + 
3Cipiip2 + Dip2, where Tpi = —Q'^ipi for z = 1, 2, upon taking the second derivative we obtain 

q^q + 2qq^ = -0^ [Ai^f + 3BiIj1^2 + 3C^i^i + DiP^] 

+ [2A'iPiiPl + 25V'?^2 + 2CV^2^i + 2DiP^il^2 + 4:BtPM2 + 4C'02V'iV'2] • 

This immediately yields 

g2 [q + n'^q] = A [(^^3 ^ 35^2^2 + SC^i^i + D^l) 

X {{Ai;f + 2BiPiij2 + CV'i)Vi + (Bipl + 2CV'iV'2 + ^V'i)V'2) 
- ((^V? + 25ViV'2 + C^p|)lP2 + {B^pl + 2C7ViV'2 + I?V'i)V'2)' 

Then after a lengthy simplification process we obtain 



■ + n\ = l [{AC - B^) (^1^2 - V'iV'2) V? + [BD - C^) (^iV'2 - V'iV'2) V; 



+ {AD - 5C)(^i^2 - M2fi^l^2\ 

= ^ [{AC - B^)i;f + {BD - C'')i,l + (^^ " BC)^^^l)2\ W\ 

where W = (V'iV'2 - V'iV'2) is the Wronskian. Let AC - B^ = BD - C"^ = -A and AD-BC = 0, 
then we obtain 

.. ^2 2AT^2 

g + f^2g = -^(VJ + Vi)- 

It is known that for harmonic oscillator Tpf + ip2 = E = const. Thus we obtain 
2AEW'^ 



q + ^^q 



The conditions on the coefficients as stated above imply that A + C = B + D = 0. 
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2.3 Mapping the Ermakov— Pinney equation to the Gambier equation 

Next we map the generalized Ermakov-Pinney equation to a reduced version of Gambler's 
equation. Let 

a 

X — 

where a is some constant parameter. By a direct computation one can show that if q satisfies 
equation p.Sp then x satisfies 



X 3^x3 - Sn^x = 0. (2.6) 

3 X 



Equation (j2.6p has a resemblance with the reduced Gambier equation. This can be easily 
generalized for time-dependent a{t). 

Proposition 3. Let x = If q satisfies the generalized Ermakov-Pinney equation (j2.5p then 
X satisfies 

4i;2 2d K o ( 2 ( a 

-:T-i-3-x3+ - 3J]2U = o. 

6 X 6 a \6\a) a J 

Proof. A direct computation of (j2.5p with q = {a{t)/x)^^^ leads to the proof. ■ 
Proposition 4. The equation 

.. A B 

q + Xq + fiq= ^ 

g-^ q 

may be mapped to the reduced Gambier equation l\1.2\) with n = —2 and a = by the following 
transformation x = a/q'^, where a is a constant. 

Proof. The proof follows a direct computation with q = ^/ax~^^'^ which gives 

^ ~ 2^ 



2 ' " 2x3/2 

Substituting these expressions into the above equation we obtain 

3x2 2A . 2B o 

X = Ax + 2/zx H ^x H X . 

2 X a 

A comparison allows us to identify the coefficients in terms of those occurring in (|1.2p with 
n = —2 and a = 0. ■ 

Note that the transformation x = a/q'^ is generally referred to as a folding transformation 
and is discussed in the next section. 



2.4 Folding transformations and the Gambier equation 

This section is devoted to the so called folding transformation of the Gambier equation. Folding 
transformations were first introduced in [22] by Okamoto, Sakai and Tsuda for the Painleve 
equations. According to their terminology folding transformations are algebraic transformations 
of the Painleve systems which give rise to a non-trivial quotient map of the space of initial 
conditions. In a more elementary sense one may say that folding transformations relate the 
solution of a given Painleve equation to the square of that of some other equation (which may 
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be the same as the initial one) [19]. One of the interesting resuhs of Okamoto et al. was the 
discovery that Piv possesses a folding transformation of degree three, mapping solutions into 
solutions of the same equation for different values of the parameters. 

Let us demonstrate the application of folding transformation for the Gambler equation. When 
6 = 0, n = 2 and a = const, the Gambler equation p.2p reduces to the following equation 

x = \-2axx X — aax . (2.7) 

2 X 2 2x ^ ^ 

Under the transformation x = —q^, this is mapped to the following equation, 

,- + 2a,2<^+{(f),+ (|)%^+(f)'l}=0, 

which is seen to be an equation of the Lienard type. On the other hand if we further set a = 
then (|2.7p reduces to 

1 ^2 q2 

X = h 2axx X , 

2 X 2 

and under the same transformation as above we find that 
q + 2aq^q+ = 0. 



,2. 

If however, we set a = but cr 7^ then (j2.7p assumes the following form 



±2 ^2 



2x 2x 

which under the transformation x = —q^ becomes 

^2 1 

4 q-^ 

2.5 Generalized Milne Pinney system and Gambler equation 

The second-order nonlinear differential equation 

q + w{tfq + Kq-^ = Q (2.8) 

describes the time evolution of an oscillator with inverse quartic potential. This is a generali- 
zation of the Milne-Pinney equation where the evolution is governed by the inverse quadratic 
potential. It is possible again to map the generalized Milne-Pinney equation to the Gambler 
equation. 

a{t) 

satisfies 



Proposition 5. Let x = If q satisfies generalized Milne-Pinney equation (12. 8p then x 



4x2 (2a n.itf\^ ^ ^ + (^V ^ 0. 



3 X V 3 a 3a2 / a^ V3Va/ a 3a^ 



Corollary 1. Let w{t)'^ = 2aa. Then the generalized Milne-Pinney equation is mapped to the 
reduced Gambier equation 
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3 Generalized Riccati differential operators 
and Gambler's equation 

In this section we construct a generalized Riccati differential operator which generates a new 
sequence of differential equations, that includes the Gambier equation. 
We begin by defining the operator 



d n — 1 X ( G b ax\ 
- -I 1 



dt n X \nx 2 n J 



(3.1) 



Here the dot represents differentiation with respect to t and n is an integer. Furthermore a and 
b are assumed to be functions of t while a is a constant. 
Next we consider the sequence of differential equations, 

G"'x = 0, m = 0,1,2,..., 

which give rise to the following ODEs: 

I a b a 2 
m = 1, Gx = =^ —X = — I — X H — X , 

n n 2 n 

which is a generalized Riccati equation. For m = 2 we have G'^x = 0, which implies 

n — 1 iP' n + 2 n — 2 x a? n 
X = h a XX + bx (T X 



n X n n X n 

„2 



+ (d — ab)x^ + ( cn a] x — ba , (3-2) 

n I nx 



where 



h 62 

'=2-T 

This procedure can be carried on leading to higher-order equations. 

In this paper we will mostly be concerned with (|3.2p . As has already been mentioned (see 
also [11]) the linearized family of equations belonging to Gambier 's classification can be obtained 
as canonical reductions of this particular equation. 

Indeed for (11.31). viz 



n— li n+2 an 9 

X = h a XX x + ax , 

n X n n 

the corresponding generalized Riccati operator is 

d n — 1 X a . 
Gi4 = X. (3-3) 

dt n X n 

This corresponds to the operator for the Gambier equation ()G14p . 

On the other hand the choice n = 1, o" = 0, a = —1 and c = 0, so that b = 6^/4, leads to the 
fifth equation of Gambler's classification scheme, namely 



-3xi; — + b[x + x^) , 



and corresponds to the following generalized Riccati operator 
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Taking the hmit n — >• oo in (|3.2p we are led to the following equation, 

• 2 

X = h axx + hx a + (a — ah)x'^ + dx — ba, 



where d = lim n{b/2 — b /A). However, it is not quite apparent what the corresponding gene- 

n— j-oo 

ralized Riccati operator for this equation is, since a direct limiting procedure on (jS.ip does not 
lead to the 13th equation of Gambler's classification. 



3.1 Hamiltonian framework for the Gambier equation (]G14p 

We consider equation (jl.3p 



n— li n + 2 a'' ^ n 

X = h a XX x + ax . 

n X n n 

When n = —2 this becomes 

X = -— + —x'^ + ax\ (3.4) 
I X I 

which may be written in the usual manner as the system 

x = y, y = - h— x'^ + ax. (3.5) 

This system can be expressed in a Hamiltonian form with the following Hamiltonian [3]: 

H{x,p,t) = ^x^ {p + b{t)f -dlogx, 

where b{t) = | /* a'^{s)ds. The canonical coordinates used are q = x andp = yx~^ — ^ J* a^{s)ds. 
The Hamiltonian equations of motion are 

■ dH 3, 

^ ^ ~dp =^ '^P^")' 

•^if 3 2/ ,x2 « 

p=-^ = -r^p^'^^-x 

Using (p + 6) = y/x'^ we find that x = y. On the other hand 

y = 3x^x(p + 6) + x^{p + b). 

Since b = l/2a^ and x = y, the last expression can be simplified to yield the second equation 
in (1331). 



3.2 Conjugate Hamiltonian equation of Gambier type 

Consider a non-autonomous Hamiltonian system in a (2 + l)-dimensional extended phase space 
{PjQjt} with a given Hamiltonian h = H{p,q,t). The equation of motion is given by 

dq_dH_ dp _ dH 
dt dp ' dt dq ' 

The key point is that the Hamiltonian must be explicitly dependent on time which yields t = 
T{p, q, h). Thus we express 

h = H{p,q,T{p,q,h)) 
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and by chain rule we obtain 



dT 


OH 


fdHy' 


dT 


dH , 


fdH\-' 


dT 


fdH 


dp 


dp ' 




dq 


dq ' 




dh 


\ dt 



Since f = ^ and 

dp _dp fdH\~^ dq _dq fdH\~^ dt _ f dHV~^ 
dh~dt\dt) ' lh~~dt\dt) ' lh~\dt) ' 

Substituting ()3.ip in ()3.3p and considering (j3.2p yields 

dp _dT dq _ dT dt _ dT 
dh dq ' dh dp' dh dh ' 

This set of equations are known as the conjugate Hamilton's equations. 
3.2.1 The Lienard equation of second type and its conjugate 

The concept of a conjugate Hamiltonian is introduced in |23] . The solution of the equation 
h = H{p,x,t), where H is a given Hamiltonian containing t explicitly, yields the function 
t = T{p,x, h). The Hamiltonian system with Hamiltonian T and independent variable h is said 
to be conjugate to the initial Hamiltonian system with Hamiltonian H. Using this construction 
Fokas and Yang [7] derived the conjugate equation for Painleve H. In particular, they presented 
a Lax pair formulation as well as a class of implicit solutions. It may be mentioned here that, 
in general, a transformation of the hodograph type does not preserve the Painleve property. We 
now apply their construction to a Lienard equation of second type. 

The Hamiltonian of the Lienard equation of second type, viz, x + f{x)iP' + g{x) = 0, may be 
written as 

H{p,x,t) = \e-'^^^\p-g{x)e'^^^h)\ 
where F{x) = J f{x) dx. Thus H = h implies 

It can be easily shown that the dual equations are 

dx^_dT^ l_^-2F{x) 

dh dp g{x) ' 

dh dx dx \g{x) ) ^ ' g{x) ^ ' 

Conjugate Hamiltonian for the reduced Gambier equation. The choice 

/(2;) = -^ and g{x) = -x^ 
leads to the conjug ate of ([331) when a{t) = ^/2 as given below 

dx _ dp 3 I— c:/9 

These equations may be solved explicitly yielding 




3 ( I h_ 

Jh + 'ciY 



(h + ci), p = =p— / ———-^dh + C2, 



where ci and C2 are arbitrary constants of integrations. 
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4 Generalized Milne— Pinney type equations 
and superintegrability 

The notion of superintegrability [6l |15] has attracted a lot of attention in recent years, due to 
its various applications and associated mathematical structures. Fris et al. [6], in 1965, studied 
Euclidean n = 2 systems which admit separability in two different coordinate systems, and 
obtained four families of superintegrable potentials with constants of motion linear or quadratic 
in the velocities (momenta). In this section we focus on a certain reduced version of the Gambler 
equation (|3.2|) which exhibits the property of superintegrability. If we set n = 2 and a = b = 
then equation (j3.2p reduces to following class of equation 



1 o ^ K 

X X + ilx = 0, 

2x X 

upon ignoring the constraint on c. Our aim is to study such an equation using a singular 
transformation of Sundman type. 

We look for nonlocal transformations, under which a given ordinary differential equation is 
linearizable. The problem was studied by Duarte et al. by considering basically a transfor- 
mation of the form 

X{T) = F{t, x), dT = G{t, x)dt. 

Here F and G are arbitrary smooth functions and it is assumed that the Jacobian J = ^^^^^ 7^ 0. 
If one knows the functional form of x{t), then the latter transformation ceases to be nonlocal. 
One of the pioneering contributors to this field was K.F. Sundman [21] who introduced the 
transformation dt = rdr in the study of the 3-body problem, where r is the dependent variable 
(radial component). In particular such transformations are especially effective for the solution 
of several nonlinear ODEs. 

Consider a second-order ordinary differential equation of the Lienard-II type 

x + f{x)±^+g{x) = 0, (4.1) 

and suppose we seek a nonlocal transformation such that it is mapped to the usual equation of 
a linear harmonic oscillator, namely, 

X + io^X = 0, (4.2) 

(here X = ^) by the nonlocal transformation 
dX 

— = A{x,t)dx + B{x,t)dt, T = t. 
X 

It is a matter of straightforward computation to show that the SODE ()4.ip is mapped to ()4.2p 
provided its coefficients satisfy the following conditions: 

+ A^- A{x)f{x) = Q, 
+ 2A{x)B{x) = 0, 

A{x)g{x) -B^ = io^. (4.3) 
Let A = Vx/v, then upon solving the first and second equations of the above set we obtain 

Vx = exp (^J f{s)ds^ , B{x) = ^. 
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The final equation ()4.3p then yields 



—9{x) - — = w . 

V 



In particular, let us set /(x) = ^, which readily yields v = ( a ^ — !)• Then A and are 
given by the expressions 

^ X = and B{x) = 

respectively while 



x4"+3 (a + 1) ■ 

Therefore, the singular Sundman type transformation is given by 



T = t, 



while the equation assumes the form 

x + -x^ + ——^^x + ^ = 0. 4.4 

X a + 1 x^'^^'^ 



Integrating ()4.2p we get 

I + w A = i (t, X, xj = const, 

where /(t, x, x) denotes the first integral. Using this recipe we obtain the first integral of (jlT 
as 

dt 



I = {a + If (^x°x + ^) + a;2 {x'^+^ exp {{a + 1)^ y 



3.2(0+1) 



As shown in [S] (and references therein) the second-order equation (14. 1|) admits a Lagrangian 
description via the existence of a Jacobi last multiplier (JLM) M(x), [13] of the form 

/a; 
f{s)ds, 

with 

1 /"^ 

L{x,x) = — M(x)x^ — l^(x), where the potential is V{x) = / M{s)g{s)ds. 

The following proposition is important in identifying superintegrable systems. 

Proposition 6. The necessary condition for mapping of V{x) to the superintegrable potential 
V{x) = AQ^ + B /Q"^ is given by the condition 

g"{x) + {f{x)g{x))' 



g'{x) + f{x)g{x)-2A 



g"'(x) + {f{x)g{x))" ^^^^^ 



g"{x) + {f{x)g{x)y ' 
where Q = y^M{s)ds and A and B are arbitrary constants. 

Clearly one can directly check that (j4.5|) is satisfied by (j4.4|) . Moreover the potential is 
mapped to 



V{x) = I x^'^gix) dx = / x2" 



cj^x (a + 1) 



31 



+ 



a + 1 x4"+3 



dx = - 



la;2(x2("+i)) l(a + l) 



2 (a + 1)2 2 x2("+i) 



Hence we have Q = x"+V(" + 1) while A = and B = -1/2. 
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5 Conclusion 

In this paper we have examined the Gambier equation and its associated systems hnking them 
to several weh- known differential equations. In particular we have established the similari- 
ty between special cases of the Gambier equation and the generalized Ermakov-Pinney and 
Milne-Pinney systems. The use of Okamoto's folding transformation leads to a reduction of the 
Gambier equation to an equation of the Lienard type. In the case of the Gambier equation ()G14p 
it is further observed that the system admits a Hamiltonian description. As a consequence we 
have studied the conjugate Hamiltonian equations for the second type of Lienard equation 
which resembles the Gambier equation in certain respects and have indicated their relevance in 
connection with superintegrability. 

Given that the extension of the Gambier system to the discrete case has been already estab- 
lished we could speculate on whether some of the results presented here could be extended to 
the discrete domain. This is a question we hope to address in some future work of ours. 
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